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THE SIMILARITY RULES FOR SECOND-ORDER SUBSONIC AND SUPERSONIC FLOW 1 

By Milton D. Van Dyke 


SUMMARY 

The similarity rules for linearized compressible flow theory 
( Qdtherl’s rule and its supersonic counterpart) are extended to 
second order. It is shown that any second-order subsonic flow 
can be related to “nearly incompressible” flow past the same 
body , which can be calculated by the Janzen-Rayleigh method. 

INTRODUCTION 

The linearized small-disturbance theory of steady com- 
pressible flow, based on the Prandtl-Glauert equation, yields 
a first approximation for thin objects moving at either sub- 
sonic or supersonic speeds. More precisely, it provides the 
first term iD an asymptotic expansion of the solution for 
small disturbances, provided that the flight Mach number 
is not too close either to unity (transonic flow) or to infinity 
(hypersonic flow). 

The similarity rule that governs linearized subsonic flow 
past general three-dimensional objects was first given cor- 
rectly by Gathert (ref. 1). It has an obvious counterpart in 
supersonic flow, and the rifles have rendered great service in 
both theoretical and experimental investigations. 

Recently, various investigators have sought to improve 
on the linearized theory by finding higher approximations 
(see, e.g., refs. 2 to 5). The second step is commonly 
referred to as the second-order small-disturbance theory, or 
simply “second-order theory.” It can be found in general 
by iterating upon the linearized solution, retaining all terms 
out of the full nonlinear equations of motion whose contri- 
bution is of the order of the square of the disturbances in 
linearized theory (ref. 3). In the simplest case of plane 
flow without stagnation points, the linearized disturbances 
are proportional to the thickness ratio r of the airfoil, so that 
second-order theory adds terms in t 2 , and higher approxi- 
mations extend the series in powers of r n . Stagnation points 
lead to the appearance of logarithmic terms, beginning with 
rHn t in the fourth approximation. The series diverges in 
the immediate vicinity of stagnation points, although it can 
be corrected there by simple techniques (ref. 6). Slender 
pointed objects, such as a smooth body of revolution, cause 
smaller disturbances than airfoils, but logarithmic terms 
always arise at the outset; hence the linearized solution 
contains terms of order r*ln r and r 2 , and the second-order 
increment then consists of terms of order r 4 ln 2 r, r*ln r, and t 4 . 
Nothing is known of the convergence of these series; they 
are perhaps only asymptotic expansions for small thickness. 
Second-order theory, like linearized theory, breaks down in 
the transonic and hypersonic ranges, though it may pene- 
trate somewhat farther into their fringes. - 

* Supersedes NAOA TN 3875, by Milton. D. Van Dyke, 1957. 


A similarity rule for second-order theory has recently 
been given in the special case of supersonic flow past thin 
flat wings by ffenain and Germain, who demonstrate its 
usefulness in theoretical studies (ref. 5). However, as in 
linearized theory, the rules for flat wings are only special 
cases of those for general three-dimensional shapes. The 
present paper is devoted to deducing the general rules for 
subsonic and supersonic flows, and examining their implica- 
tions. In particular, it is shown how the rule for subsonic 
flow relates the second-order solution for any object to nearly 
incompressible flow past the same body, which can be calcu- 
lated by the Janzen-Rayleigh method. 

The author is indebted to Wallace D. Hayes for suggesting 
several improvements that have been incorporated in the 
present version of this paper. In particular, the procedure 
for recovering the second-order solution from the Janzen- 
Rayleigh solution (p. 930) is simpler and more logical than 
that originally given in NACA TN 3875. 

DERIVATION OF RULES FOR BODIES OF REVOLUTION 

A body of revolution is the simplest shape that is not a 
special case, but displays the full generality of the existing 
s imil arity rules for subsonic, supersonic, transonic, and hy- 
personic flows. The same can be shown to be true of the 
second-order rules to be discussed here. Hence for clarity 
of exposition, the second-order rules will be derived in 
detail only for an axisymmetric body at zero angle of 
attack. The rules for general three-dimensional thin or 
slender objects will thereafter be stated without proof. The 
subsonic and supersonic cases are so similar that they can 
be treated simultaneously. 

Let the body be described by r—rR{x), where r is a thick- 
ness parameter or characteristic slope (say, the maximum 
slope, average slope, thickness ratio, or the like), and R(x) 
is a function of order unity (fig. 1). As usual in similarity 
analysis, the characteristic slope r is regarded as a parameter, 
so that different values of r correspond to affinely related 
members of the same family of bodies. 
r 

| r — r*rR(x) 



Figure 1 . — Notation for body of revolution. 
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To second order the flow is irrotational, so that there 
exists a velocity potential $(x, r; M, y, r). This notation 
indicates that for each family of bodies (associated with a 
given function R(x)), the flow field is regarded as depending 
not only upon the two independent variables x and r but 
also upon the three parameters following the semicolon: 

M free-stream Mach number 
y adiabatic exponent of gas 2 

r characteristic slope of body 

The aim of a similarity analysis is to transform the prob- 
lem so as to reduce the number of parameters appearing in 
it. If that can he accomplished, flows having different 
values of the original parameters are related provided only 
that the reduced parameters are equal. The transforma- 
tion to be used here consists in separating the dependent 
variable $ into several components, and then stretching 
each component and the independent variables by factors 
that depend upon the original parameters. It is convenient, 
and involves no loss of generality, to leave streamwise 
coordinates unchanged, so that r is to be stretched but not x. 

Perturbation potentials are first introduced by setting 

■fr-—x-\-^ J r(p-\- ... (1) 

V-/ CD 

where <t> is the potential of linearized theory, and <p the 
second-order increment. 

RULES FOR LINEARIZED THEORY 

The linearized problem is 

□ ^= (1— HI 1 ) <l>xz+<l>rr+~-—0 

<£->0 at infinity ” ^ 

<i> r =rR'(x) at t=tR(x) „ 


from two to one is (temporarily suppressing the dependence 
on parameters) 

<t>{x,r)=^F{x,p) (3a) 

and 

P=0r (3b) 

where 

{ Vl —M 2 for subsonic flow 

, (3c) 

Vilfi-l for supersonic flow 

Then the problem becomes 

AF=F„+^±F a =0 

p 

F^O at infinity ^ 

F p =PtR'(x) at p=0tR(x) , 

where here and later the upper and lower signs apply, 
respectively, to the subsonic and supersonic problems. 

The transformations of <t> and r have been so chosen that 
the problem is reduced to one involving the two parameters 
M and r not separately, but only in the combination /Sr. 
This is the similarity parameter. Two subsonic or super- 
sonic flows past bodies of the same family are related if 
the corresponding Mach numbers are such that the param- 
eter fir is the same for both. The nature of the relationship 
is found by reintroducing the dependence on parameters 
into equation (3a), which gives the similarity rules 

<t>(x,r;M, t ) F(x, /3r; /Sr) (5) 

SECOND-ORDER RULES 

Thee second-order problem is found to be (ref. 3) 


The first relation is the linearized Prandtl-Glauert equation. 
The second is a statement, sufficiently definite for present 
purposes, of the requirement that the flow approach a uni- 
form stream far from the body in almost all directions. The 
third is the linearized condition of tangent flow at the body 
surface. 3 The linearized problem is seen not to involve y, 
so that the solution depends upon only the two parameters 
M and r. 

The similarity rules can be obtained by considering arbi- 
trary scale transformations of 4> and r. It is readily found 
that, the only choice that reduces the number of parameters 


~0<f>=M i { [(y+l)M t -\-2(l—M*)]<l> I 4> xx -\-2<f> r <}> ir -{-4> r i <l > rr } ' 
<p—*0 at infinity 

tp r —T<f> x R'(x) at r=rR{ x) „ 


( 6 ) 


Note that the first equation contains not only quadratic 
terms on the right-hand side, but also the triple product 
4> T 2 4>rr whose contribution is of second order in some cases. 

The parameter y appears only linearly in the combination 
( 7 + 1 ) and can accordingly be separated out. Thus tlio 
appropriate transformation is found to be 


* Hayes has pointed ont (ref. 7) that to second order an imperfect gas corresponds to a 
polytropic gas having a y equal to the free-stream value of 

m— + m. 

where c Is the speed of sound and p the density, the partial derivative being taken at constant 
entropy «. 

* In what Is generally called the slender-body approximation, the body Is assumed to be 
so smooth and slender that the tangency condition can be Imposed on the axis rather than 
on the actual surface, for bodies of revolution in the form lim T<h=r*R(z)R r (jc). Thus 

7—+0 

slender-body theory is a further approximation within linearized theory (being. In fact, the 
leading term in the asymptotic expansion of the linearized solution for small thickness r). 
Consequently, the slender-body solution obeys the similarity rules of linearized theory, and 
he second-order slender-body solution likewise obeys the second-order similarity rules. 


v(x,r)=± \j i (x >P )+M 2 f i (x, P )+(y+l) ^ / 3 (a*, P )] (7) 


Then equating like powers of M 2 yields the following set of 
three problems for / 1 / 3/3 in which the parameters M, 7 ,t 
appear again only in the form of the single similarity param- 
eter /St: 

A/i=0 "I 


fi — >0 at infinity >- 

f ip =fiTF x R'(x) at p=firR(x) J 


(8a) 
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Af i =±2F x F xx +2F p F xp +FfF„ ' 

/ 2 -> 0 at infinity >- (8b) 

f 3p = 0 at p=firR(x) 

A/, =F x F a 

f a -> 0 at infinity >- (8c) 

^=0 at p = firR (x) , 

Then reintroducing the explicit dependence on parameters 
into the functions of equation (7) gives the similarity 

rules for the second-order increment in perturbation po- 
tential: 

<p(x,r,M, [/!<*,<&•; fir)+M%( ) + (t+1) ^ M )] (9) 

where the arguments of / 2 and f 3 are the same as those of 
fi. Hereafter the arguments of later terms will be omitted 
in this fashion when they are the same as for the leading 
term. It should be emphasized that, just as in linearized 
theory, the subsomc and supersonic rules are quite distinct, 
although they have the same form (9). Because of the 
different definitions of fi, and the resulting ± signs in equa- 
tions (4) and (8b), a supersonic flow is not related to a sub- 
sonic flow. Discussion of these results is deferred to the 
general case. 

RULES FOR GENERAL BODIES 

Consider a family of general three-dimensional bodies, 
whose members aro derived from' one another by a uniform 
magnification or reduction of all dimensions normal to the 
free stream (fig. 2). Each member of such a family can be 
characterized, as before, by some characteristic slope r. It 
may be emphasized that r can be identified with thickness, 
camber, or angle of attack, all of which vary together for 
related bodies. 



Figure 2. — Example of two related bodies. 


The preceding analysis can be extended in a straight- 
forward way to such general bodies, at the expense only of 
typographic complexity. Both cross-stream dimensions 
behave in the way that r did before. Hence the subsonic 
and supersonic second-order rules for the velocity potential, 
corresponding to equations (5) and (9), are, in Cartesian 
coordinates: 

jj- 4- (x,y 1 z-,M,y,T)=x-\~ F(x,fiy,fiz; fir) + 

± [j f 1 (x,fiy,fiz;fiT)+M*M ) + (t+D jpU )] (10a) 

Differentiation yields the corresponding rules for velocity 
components (those for w having the same form as for v) : 

jj-=l-\-^U( x ,fiy,fi z ',fiT)-\- 

^u^fiy^z-fi^+MHiai ) + (y+l)~^th( )J (10b) 

v ( x >fiv, fi z ]fi T ) + 

jp [j>! (x, fiy, fiz;fir) +M % ( ) + (y+l)^®s( )J (10c) 

(The functions appearing here are actually related to deriva- 
tives of the functions in equation (10a), but the connection 
is of little interest.) To second order the pressure coefficient 
is given by 

O p = —2 <p x — (<j> „*+&“) —2<p x —2(4>,<p p J r <t>z<Pz) + (M 5 — 1) ^**4- 

where the terms in the second line may be significant for 
slender shapes. Substituting the expressions (10) for veloc- 
ity components and simplifying shows that the similitude for 
pressure coefficient has the same form as that for the stream- 
wise velocity increment Au/U a : 

O p (x,y,z;M,y,T)=^P(x,fiy,fiz;fiT)+^p l (x,fiy, fiz; fir)+ 

7+1) ^o] (10d) 

The similarity rules for the perturbation stream function 
in plane flow are the same as those for v/U w (eq. (10c)). 

ALTERNATIVE FORMS 

As with other similarity rules, an unlimited number of 
alternative forms can be produced by multiplying by powers 
of the similarity parameter. Thus, of the many possible 
alternatives to the second-order rules (lOd) for pressure 
coefficient, two of the most useful are: 

C P =^F(x,fiy, fizifi^+^jpd )+M>p 2 ( )+(y+l) ~p 3 ( )] 

(lOe) 

C p =T'P(x,fiy,fiz;fiT)+T*[M)+M>p a ()+(v+ l)^ft( )] 

(10f) 
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In addition, the first two second-order terms can be manip- 
ulated, using the connection between M 1 and /S 2 , to yield 
additional alternative forms such as the following, which 
correspond to the three forms above: 

0 P =^P(x,Py,pz;M+^i[pi( )+^rJp2 f ()+(y+l)^~rPz( )] 

(10 g) 

(lOh) 

C,=t*P(x, p y,fiz; /3r) +r< [/3*&( ) +W ( ) + (r+1) U )] 

(lOi) 

FORCE COEFFICIENTS 

The rules for pressure imply rules for the lift and drag 
coefficients. For example, equation (lOe) leads to 

0 L (M, y,r) = T - L(dr) 4~ [U (dr) +M*h(P r) + fy+1) ~ *»(j8r)] 

(10j) 

G D (M,y, r) P(dr)+^ [d t (dr)+M^(dr) +(y+l) ^d^r)] 

(10k) 

if the coefficients are referred to some plan-form area. If 
some cross-sectional area is used, each term is reduced by one 
power of r. Various alternative forms are again useful. In 
the case of lift coefficient, one will ordinarily choose to 
identify r with the angle of attack. 


RULES FOR QUASI-CYLENDRIGAL BODIES 


A special class of objects must be distinguished, which will 
be termed quasi-cylindrical bbdies. These are shapes that 
lie everywhere so close to some cylinder (not necessarily 
circular) parallel to the free stream that to a first approxi- 
mation the condition of tangent flow can be imposed at the 
cylinder rather than on the actual body surface. likewise, 
in second-order theory the tangency condition can be trans- 
ferred to the cylinder by Taylor series expansion. The sim- 
plest example is an airfoil whose thickness, camber, and angle 
of attack are so small that the tangency condition can be 
transferred from the airfoil surface to a mean plane parallel 
to the stream (fig. 3). Another example is an open-nosed 



/ / / / 7 /7 / 

Figure 3. — Examples of quasi-oylindrieal bodies. 


body of revolution whose radius varies only slightly. Others 
are biplanes, cruciform wings, any of these in an open or 
closed wind tunnel, in combination with one another, etc. 

A quasi-cylindrical body can be regarded as consisting of 
a skeleton upon which is superimposed a small slope distri- 
bution. The skeleton is simply the projection of the body 
onto the basic cylinder. For example, the skeleton of the 
quasi-cylindrical body of revolution is the circular tube 
shown dashed in figure 3. ' 

The special place of quasi-cylindrical bodies in similarity 
theory arises from the fact that the skeleton and the slope 
distribution can be varied independently. This extra 
freedom is important. For example, it leads to a useful 
transonic similarity rule for quasi-cylindrical bodies whereas 
none exists for general shapes. It is convenient always to 
leave streamwise dimensions unaltered. Hence, we consider 
families of quasi-cylindrical bodies that are derived from 
one another by a lateral compression or expansion of the 
skeleton, and a quite independent magnification or reduction 
of all surface slopes. Two members of such a family are 
shown in figure 4. 




Figure 4. — Example of two related quasi-cylindrical bodies. 


Distortions of the skeleton will be measured by some 
characteristic “aspect ratio” A. It is important to note 
that the term “aspect ratio” is used here in a very general 
sense to mean any typical ratio of gross cross stream to 
streamwise dimensions. For example, in the last shape 
in figure 3, the ratio of wind-tunnel height to airfoil chord is 
an appropriate characteristic aspect ratio. Changes of 
slope are measured, as before, by some characteristic slope r. 

The preceding similarity rules can be simplified for quasi- 
cylindrical bodies by using the facts that first-order per- 
turbation quantities are directly proportional 4 to r, and 
second-order terms to t 2 . The simplification can be earned 
out by first imagining the quasi-cylindrical body to be 
restricted to be a general body, which means that both 
pA and /Sr must be the same for similarity. Then consider 
the preceding rules in the particular alternative forms in 
which r appears explicitly outside the first-order term and 

* This ts by no means true lor general bodies; as noted previously, the llrst-order pressure 
coefficient on a smooth slender pointed body of revolution varies as r*/n t for small r. 
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r 2 outside the second-order terms. For the pressure co- 
efficient, this form is that of equation (lOe) : 

0„=lp(xMPz; P tM)+ : p [ft( )+m*m )+(y+i) )] 

So far the functions P,Pi,P 2 ,P 3 have been supposed to depend 
parametrically upon both PA and /Sr. However, the first- 
and second-order terms can be proportional to r and P, 
respectively, only if the supposed dependence upon /Sr is 
nonexistent. Hence, the similarity parameter is /S A alone, 
and the rules for pressure are (dropping bars from the 
functional symbols) : 

O p (x,y,z; M, V,r,A)=^P (x, fry, pz; pA) ^p l (x,Py,Pz;pA)+ 

M%()+(y+ l)^jp,()] (11a) 

The coiresponding rules for the potential and velocity 
components can, if desired, be written down by inspection 
from equations (10). 

With p 2 =pi±pi (where, as before, the upper sign applies 
to subsonic and the lower to supersonic flow) , these rules can 
be rewritten as 


0 P =~P( )+P [*,( )+(y+ 1) ~P3 { )] (lib) 

and this is the result that Fenain and Germain found in then- 
treatment of the flat diamond cone in supersonic flow (ref. 5). 


CONNECTION WITH HAYES’ RULE 

For plane flow past a single body, Hayes has discovered a 
remarkable rule for the second-order surface pressure (ref. 7). 
It implies that, on the surface, the functions in equation 
(11a) are such that pt—0 and pi=ip 3 . Hence, 


q =IjP(x)+t 3 . fr+l)M 4 -H( j p ( x \ 

W. 4(1— M 2 ) 2 Vl( ) 


( 12 ) 


In supersonic flow this is simply Busem arm’s well-known 
second-order solution, P being twice the local slope of the 
surface and Pi twice its square; in subsonic flow P and p x are 
more complicated (ref. 8) . This rule implies a corresponding, 
but more complicated, rule for surface velocity (ref. 8) . 

In addition to the restriction to single bodies and plane 
flow, these rules are not s imil arity rules in the sense of the 
preceding results, because they apply only at the surface 
rather than throughout the field. 

EXAMPLES 


The rules will be illustrated by two simple examples, at- 
tention being confined to the surface pressure coefficient. 

SLENDER CIRCULAR CONE IN SUPERSONIC FLOW 

Broderick has derived the second-order slender-body 
solution for a circular cone at zero angle in a supersonic 


stream (ref. 2). The surface pressure coefficient on a cone 
of slope t is 

o,=p(2;»!-i)+^(;»J : )‘- 

(5M>-1) is^+y^'+i+w+lj'lr] (13) 
This has the form of equation (lOi) with 

t>= 31»*|:-5i»|+f ?,-l 

WAVY WALL IN CLOSED SUBSONIC WIND TUNNEL 

Consider the sinusoidal wall y=r sin a; at a distance h from 
a flat wall (or a distance 2 h from its mirror image) as indi- 
cated in figure 5. Subsonic flow between the walls at a mean 



/■ 




\ 


\ ^ 





\ 


y 

/ n / / / / - /, L — L — L — L L L — L L 



Mach number M can be readily calculated to second order 
by separation of variables. The resulting pressure coeffi- 
cient on the surface of the wavy wall is 


< 7 , 


2t l+e-** . . , (y+1)M‘+4(1-M 2 )^ 

t=«s sm z+r 2 ^ X 


”• p 1 — e-W 


4(1 -M 2 ) 2 

l+le-w+e-w 


M K 

2 ^ +1 )^ t5 (I +7 


(l-e- 5 ^) 2 

phe-w 


cos 2z— 2r 2 




(l-e-wy 


(l-e-wy * 

cos 2x (14) 


The relevant aspect ratio is the height h (which is really a 
multiple of the height-chord ratio, because of the choice of 
scale for the wavy wall) . Thus the result is seen to have the 
s imili tude of equation (11a). As the tunnel height increases 
indefinitely, the last two terms disappear, and the remainder 
follows the similitude of equation (12) for the surface of a 
single plane body. 
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REDUCTION OF SUBSONIC PROBLEM TO NEARLY INCOM- 
PRESSIBLE FLOW 

In linearized theory, an important application of the 
similitude is Gothert’s rule, which reduces any subsonic flow 
problem to a related incompressible flow (ref. 1). As the 
rule is usually stated, the incompressible flow is that past a 
thinner affinely related body. However, the incompressible 
solution for one member of an affinely related family of bodies 
determines that for all other members, so that the subsonic 
flow may, if desired, be related to the same body rather than 
a thinner one, and that viewpoint will be adopted here as 
being the simplest. 

In second-order theory, the explicit appearance of terms in 
M 3 and (y+1)A/ 4 in equations (10) means that reduction to 
an incompressible problem is impossible (except for the 
special case of the surface of a single plane body, where eq. 
(12) applies). The second-order problem can, however, be 
reduced to a nearly incompressible flow. 

Flows at low Mach numbers can be calculated by the 
Janzen-Rayleigh method, which involves iterating on the 
incompressible solution to obtain a power series in M ?. 
Thus the velocity potential is approximated by 

~ $ (x,y, z;M,y, r)=$ 0 (*,y, 2; r) +MH> l (x,y, z; r) + 

(7+1) )+ . . . (15a) 


the first-order term F and the second-order increment j\\ 
combining them as F t =F-\-fi gives 

)+(T+l) I? /.( ) (10m) 


The Janzen-Rayleigh solution is now to be manipulated 
into this same form. The first three terms as given in equa- 
tion (15a) are equivalent to 


1 M i 

jj- H^,y,z;M ) y,T)=^ 0 (x,y ) z;r)+-^- ) + 

(T+l)f *,()+o(f) (16b) 

which may be rewritten as 

Q—UcX , . . M 2 , . . , . M 4 . s 1 „ fM A \ 

- <pi{ ) + (Y+l) -p <Pi( )+0 \JjTJ 

(16c) 


with 


<Pa= 


^=- T 2 


<P*—^ 


_$o— x') 


$1 


(ISd) 


The two terms in M* are ordinarily considered together, 
but for present purposes it is essential to separate them 
because only is required. This is fortunate because 
can be calculated almost as easily as $j, whereas the deter- 
mination of $ 3 is much more difficult. 

The small-disturbance and Janzen-Rayleigh series repre- 
sent two different asymptotic expansions of the actual solu- 
tion. They are believed to complement each other, so that 
an expansion of the Janzen-Rayleigh solution for small 
thickness must be identical with the expansion of the small- 
disturbance solution in powers of M s , as has been verified 
in all worked examples. This fact permits the small-dis- 
turbance solution to be recovered from the Janzen-Rayleigh 
series. The converse is not true, however, except for bodies 
without stagnation points, because the small-disburbance 
expansion is not uniformly valid near a stagnation point. 

PROCEDURE FOR RECOVERING SECOND-ORDER SOLUTION 

Another alternative form for the second-order velocity po- 
tential of equation (10a), which is useful here, is 

jj- $ (x,y,z;M,y, r) =x+r 3 F(x,Py, £ 2 ; /Sr) + 

V s [/i( )+^/.( )+Cr+i) f, 3 ( )] 

Xioo 

(Here F and f u f 2 , f 3 are not the same functions as in equa- 
tion (10a) but related ones; in the notation of equations 
(lOh) and (lOi) they are actually F and/i, / 2 , / 8 .) 

For presen t purposes it is unnecessary to distinguish between 


Finally, with the aid of 

2 /=(¥y/i+^=^[ 1+ I yr+° (fr)] ( 16 ) 

and corresponding expansions for 2 and r, this may be 
re-expressed as 


=W.(»,/gy,gg ; j3r)+^[ Pl ( )+f 1 ) + 
^? 5 o |9j ()+^-% r ()J+(Y+l)-^-^() (15e) 


which is the desired form. Here, for example, ^(r,/S?/,/S 2 ;/Sr) 
means (d/dy)po(a;,y, 2 ; r) evaluated at x=x, y=Py, z=/3s, 
and r=/Sr. 

The second-order solution is thus recovered from the 
Janzen-Rayleigh approximation simply by calculating in 
turn the expressions in equations (15d) and (15e). The 
procedure caD actually be expressed by a single equation 
as follows. From the Janzen-Rayleigh approximation in the 
form of equation (15a), the second-order small-disturbance 
solution is recovered according to 


$ 1 1 f $0 {x,y,z-,T)—x , Af 3 r$i 0 -x \ , 

E7. + \ r s r* A 




z-)pz 

r-*j9r 


(17) 


where i> n means § n {x,y , z-,t) throughout, and subscripts 
indicate differentiation. 
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APPLICATION TO PARABOLA 

As an example, consider plane subsonic flow at zero angle 
of attack past the parabola described by y=r-yfx (fig. 6). 
y 



Figure 0. — Parabola at zero angle of attack. 


The Janzen-Rayleigh solution including terms in (y+1 )M i 
has been calculated by Imai (ref. 9). Although the velocity 
potential is complicated, it simplifies when only second- 
order terms in r are retained to 

|AZV ln^~— 


Mv [ 4 (T+W 1+, “^ + 


Here £,17 are parabolic coordinates related to the Cartesian 
coordinates by 


so that to second order 


2 

2ij 2 =V x 1 -jry i —x J 

In this case the expressions given by equation (15d) are 
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This example illustrates the fact that for planar systems 
these terms are not of order unity in r. Then according to 
equation (15e) the second-order small-disturbance solution is 

$ ,r. 1M S ,/, ? , ?+?\ 

U~ x+ ^ 8 t? T AV j 

7+1 ihf 4 ,r rj 4 . j £ 2 +^ s 1 , n . 

~ W ~ p~ r 1 (?++) 1 1 lTb ~ W ^ r \ (19a) 


where 


2?=^+py+x 

2 5p=y* 2 +jsy-z 


This result is of interest because it apparently cannot be 
found directly. Plane small-disturbance flows can be calcu- 
lated easily if one adopts the thin-airfoil approximation of 
transferring the boundary conditions to the line y=0 by 
Taylor series expansion, but that approximation fails near 
round noses in the second approximation and, as a con- 
sequence, divergent integrals arise (ref. 8). Instead, one 
can try to treat the round nose more carefully using con- 
formal mapping (cf. ref. 11, pp. 361-367), but the result is 
found to be indeterminate to the extent of a multiple of 
ln(£ 2 +7\)- This is the potential of a point source at the 
origin, which is an eigensolution, the proper multiple of 
which (appearing in eq. (19a)) is not determined by the 
suggested method. 

The second-order increment in equation (19a) is seen 
to include terms in T^lnr, whose function is to render the 
argument of the logarithm dimensionless. However, these 
terms are simply constants, so that no logarithms of thick- 
ness appear in the actual flow quantities such as velocity 
and pressure. As remarked in the introduction, logarithmic 
terms in thickness arise in the actual flow disturbances only 
in the fourth approximation. 

The second-order small-disturbance solution for the stream 
function can in the same way be extracted from Imai’s 
Janzen-Rayleigh solution, and the result is found to agree 
with that calculated directly by Kaplan (ref. 10) using con- 
formal mapping. It contains no terms in In r. (The direct 
approach succeeds for the stream function, although it fails 
for the velocity potential, because the tangency condition 
is imposed on the mass flux, which is affected by the above 
eigensolution.) Then using the connections between the 
stream function and velocity potential, one can verify the 
correctness of equations (19). 

CONCLUDING REMARKS 

UTILITY OF THE RULES 

The second-order rules are scarcely suited for correlating 
experimental clata, since tests on four related bodies would 
be needed in order to isolate the four functions involved. 
That they are, however, useful in theoretical analyses has 
already been pointed out by Fenain and Germain in the 
special case of supersonic flow past flat wings (ref. 6). Pre- 
vious investigators had calculated (erroneously, as it turns 
out) the second-order solution for the flat diamond cone 
shown in figure 7, and carried out numerical computations 



■v. ^ 

Figure 7. — Flat diamond cone. 



932 


REPORT 1374 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


for three different Mach numbers and four values of ,the 
parameter tan X (reported in ref. 11). Because the latter 
is the similarity parameter jS A of equations (11), failure to 
take advantage of the similitude resulted in three fold un- 
necessary duplication of computing labor. 

The reduction to nearly incompressible flow assumes im- 
portance for bodies with stagnation points. The small- 
disturbance assumption is violated, and, as was noted in 
the example of the parabola, the second-order solution 
consequently cannot be found directly. For bodies of 
revolution the difficulties appear to be even more severe. 
In such cases it is convenient to calculate the Janzen-Ray- 
leigh solution, and from it extract the true second-order 
solution by the procedure outlined above. This process 
has been carried out for the paraboloid of revolution in 
reference 12. 

NONUNIFORMITY IN SUPERSONIC FLOW 

It should be noted that in supersonic flow the s imil arity 
rules may fail in localized regions. Shock waves or expan- 
sion fans spring from comers and edges, and in their vicinity 
the formal interation procedure is not uniformly valid. The 
similarity rules for surface pressure fail locally when such 
waves intersect other parts of the body, as in the case of a 
triangular wing with leading edges ahead of the Mach cone. 
The rules for integrated lift and drag are correct to first 
order, but may be in error in second-order terms. These 
difficulties can in principle be eliminated by straining the 
coordinates according to LighthilTs technique (ref. 13) . 

FURTHER EXTENSIONS 

The similarity rules can readily be extended to third and 
higher order in the same fashion (except for additional com- 
plications in supersonic flow because of the ultimate appear- 
ance of significant vorticity engendered by curved shock 
waves). The similarity parameter remains unchanged; the 
complexity arising in a proliferation of functions multiplied 
by powers of ( 7 +l)’“ilf 2 "/S _J '’. Likewise, the small-disturb- 
ance solution to any order can be recovered from the nearly 


incompressible solution provided by an appropriate number 
of terms of the Janzen-Rayleigh solution. 

Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Oct. 18, 1956 
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